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THREE-FREQUENCY NUCLEAR QUADRUPOLE RESONANCE (NQR) 



Technical Field 

The present invention relates to nuclear quadrupole resonance (NQR). More specifically, the present invention 
relates to the excitation of two frequencies to create the desired signal at a third frequency, and also the method and 
techniques for the creation of spin echoes and other signals analogous to those obtained with standard multipulse 
single-frequency NQR methods. 

Background Art 

Nuclear quadrupole resonance (NQR) is a technique for detecting target specimens containing sub-kilogram 
quantities of narcotics and explosives having quadrupolar nuclei. Such substances include nitrogenous or chlorine- 
containing explosives and narcotics. Basically, quadrupolar nuclei will exhibit nuclear quadrupole resonance-a 
change in the angle of nuclear spin with respect to its quantization axis when it is excited by radio frequency (RF) 
radiation pulses at a particular frequency. In the better-known nuclear magnetic resonance (NMR), the quantization 
axis is determined externally by the direction of the applied magnetic field. In NQR, the quantization axis is 
determined by molecular parameters. As with NMR, different chemicals require pulses of different frequencies 
(i.e., different nuclear quadrupole resonance frequencies) to cause precession in nuclei. A device used to detect 
magnetic or NQR resonance in the quadrupolar nuclei of a target specimen is tuned to emit pulses at the frequency 
corresponding to the resonance frequency of the nuclei desired to be detected. A typical NQR excitation/detection 
circuit consists of an inductor capacitively tuned to the NQR frequency and nominally matched to the impedance of 
a transmitter or receiver by another capacitor or inductor. In this regard, the present invention is related to other 
methods of NQR detection as taught in U.S. Pat. No. 5,206,592 issued Apr. 27, 1993 to Buess et al. for 
DETECTION OF EXPLOSIVES BY NUCLEAR QUADRUPOLE RESONANCE, and U.S. Pat. No. 5,233,300 
issued Aug. 3, 1993 to Buess et al. for DETECTION OF EXPLOSIVE AND NARCOTICS BY LOW POWER 
LARGE SAMPLE VOLUME NUCLEAR QUADRUPOLE RESONANCE, and U.S. Pat. No. 5,365,171 issued 
November 15, 1994 to Buess et al. for REMOVING THE EFFECTS OF ACOUSTIC RINGING AND 
REDUCING THE TEMPERATURE EFFECTS IN THE DETECTION OF EXPLOSIVES BY NQR, and U.S. 
Pat. No. 5,608,321 issued March 4, 1997 to Garroway et al. for A MEANS FOR DETECTING EXPLOSIVES 
AND NARCOTICS BY STOCHASTIC NUCLEAR QUADRUPOLE RESONANCE (NQR), and U.S. Pat. No. 
5,804,967 issued Sept. 8, 1998 to Miller et al. for A MEANS FOR GENERATING SHORT RF PULSES WITH 
RAPID DETECTOR RECOVERY IN STOCHASTIC MAGNETIC RESONANCE, and U.S. Pat No. 6,242,918 
issued June 5, 2001 to Hepp et al. for APPARATUS AND METHOD FOR REDUCING THE RECOVERY 
PERIOD OF A PROBE IN PULSED NUCLEAR QUADRUPOLE RESONANCE AND NUCLEAR 
MAGNETIC RESONANCE DETECTION SYSTEMS BY VARYING THE IMPEDANCE OF A LOAD TO 
REDUCE TOTAL Q FACTOR, and U.S. Pat. No. 6,054,856 issued April, 25 2000 to Suits et al. for MAGNETIC 
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RESONANCE DETECTION COIL THAT IS IMMUNE TO ENVIRONMENTAL NOISE, and U.S. Pat. No. 
6,104,190 issued August 15, 2000 to Buess et al. for MEANS FOR DETECTING NTTRAMINE EXPLOSIVES 
BY ,4 N NQR OF NTTRO GROUPS. 

, ., It is common to detect a magnetic resonance signal by placing a sample to be measured in a tuned, electronically 
resonant tank circuit. Then, the response of the tank circuit to the electromotive force produced by nuclear or 
electronic spins in the sample is measured. With Nuclear Magnetic Resonance (NMR) or Nuclear Quadrupole 
Resonance (NQR), the sample is placed in or near an inductor, commonly referred to as a coil, that detects AC 
magnetic fields. The inductance of the coil is tuned with a parallel and/or series capacitance to make the circuit 
electrically resonant at the measurement frequency. One or more additional reactive impedances (inductors or 
capacitors) are typically added to adjust the resistive impedance at resonance to a particular value which optimizes 
the detection sensitivity. 

Although the NQR detection technique works reasonably well in some circumstances, one of the challenges that 
NQR faces is that when applying the RF magnetic field needed to detect NQR, one can create false signals at or 
near the NQR frequency. For example, acoustic vibrations are created in certain magnetic metals. In turn as the 
magnetic domains in the metal vibrate, they induce a signal back in the NQR receiver coil at essentially the same 
frequency as the driving frequency. For example, if one has a suitcase with chrome trimmings, the detected signal 
may indicate that there are explosives in the suitcase when in fact it is acoustic ringing from the chrome that is 
observed. 

A multiple-frequency technique would eliminate the false alarms due to acoustic ringing. For example, for a spin-1 
nucleus such as "N, the three transition frequencies between the levels are discrete. The basic concept is to use 
two of the frequencies to excite the third transition, and to then detect the signal from this third transition. This 
avoids any interference from acoustic ringing since the RF is not applied at the frequency that is detected The 
method disclosed here is the first direct observation of a NQR transition near a frequency that has not been used for 
irradiation. 

There is another advantage to not irradiating at the observation frequency: the probe recovery time is greatly 
reduced. Following a high power RF pulse one has to wait for the energy stored in the coil to dissipate before 
signal can be detected. This recovery time is typically about 20 coil ringdown time constants. The ringdown time 
is proportional to the coil quality factor, Q, and inversely proportional to the frequency. At low frequency a 
substantial portion of the signal may be lost during the recovery time. In three-frequency NQR, the energy stored 
in the observation coil is limited to that which may leak in from the excitation at other frequencies. Therefore 
improved sensitivity is expected with three-frequency NQR at low frequencies when the signal lifetime is short 
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The nuclear wave function evolves under a Hamiltonian consisting of the large time-independent quadrupole term 
and the much smaller time-dependent terms corresponding to the alternating RF magnetic fields applied at two of 
the three characteristic NQR frequencies. For example, with frequencies for spin-1 nuclei and the principal axes 
frame (x, y, z) of the electric field gradient tensor at the quadrupolar nucleus, the quadrupolar Hamiltonian Hq is 
H Q = e 2 qQ[(3ll - / 2 ) + Vlil] - /*)]/4 where / is the nuclear angular momentum operator, ^the 
asymmetry parameter of the electric field gradient, q the field gradient, and Q the quadrupole moment of the 
nucleus. The transition frequencies between the eigenfunctions of H Q are CO ± = (3 ± rj)e 2 qQ /A h 

and COq = T]e 2 q Q /l fi . Here, it is assumed 77 * 0 or 1 to avoid degenerate energy levels and transition 
frequencies which complicate the calculation. To be concise, only the observation at a> + arising from irradiation at 
co. and GJb is here treated in detail, but similar results are expected for other three-frequency combinations. 

For a spin-1 nucleus each transition is allowed under a different orientation of the applied field in the principal axes 
frame (i.e., (+|/|-> = <+|/J->, <0|/|-> = <0|/ y |->, and <+|/|0> = <+|/,|0». Therefore, for a single crystal, the most 
efficient NQR excitation and detection occur when an RF magnetic field of frequency g% is applied along the z- 
direction, that of a> m along y, and the 0+ detection coils are sensitive to magnetization oscillating along x. Similarly, 
for a powder sample it can be shown that the maximum signal is obtained if the two RF magnetic fields are 
mutually perpendicular to one another in the laboratory frame (x\ y\ z'). The received 'three-frequency' NQR 
signal then arises from a magnetization which is orthogonal to both the applied RF fields. 

The Hamiltonian for the interaction of the nucleus with an RF pulse of magnetic field strength B lm and frequency co. 
applied along the x'-axis is H u = -hy u B x „ / v cos cot and with an RF pulse of strength B l0 and frequency tf^ along / 
is H lQ = -hy N B l0 / y .cos cotf. ( is the gyromagnetic ratio of the nucleus.) The lab frame operators can be expressed 



as 



/ T , = (cosacosjScosy- sin«siny)/ x + (sin acos/3cosy+cosasin y)l y - sin£cosy/ 2 
/ = (-cosacos^siny- sinacosy)/ x + (- sin acos^sin y + cos acosy)l y + siny3siny/ z 

where a, p, and y are Euler angles describing the relative orientation of the principal axes and lab frames. If it is 
assumed each RF pulse only excites one transition, H 10 and H^can be simplified as 

H x _ = -h Yn b i- ( sin acos/3cosy+cosasiny)/ y coscoj = -nQ { _ l y coscoj 

(2) 

H l0 =-n y N B 10 (sin £sin y) 1 1 cos co 0 t = - hQ l0 / z cosco Q t 



where the dependence on crystal orientation is now contained implicitly in the newly defined terms Q,„ and Q,_, 
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the effective RF-nutation rates. 

Using the above Hamiltonians, the wave function 1 5Fi»> is found after a single RF pulse of duration r p at a 
frequency can be written as a simple rotation of the original wave function around z, 

, where 0= Q 10 1„ 12. Similarly an RF pulse at the frequency co_ is equivalent to a 
rotation about y by Q,. r„ 12. Furthermore, the wave function after simultaneous irradiation at u and H can be 
shown to be equivalent to a rotation about an axis in the y-z plane, rotated from the z-axis by an angle £ where tan 

{ . o. Q , ™ i*w> = m«» . where vsTO77,, /2 lste 

angle of rotation defined by the effective RF-field generated by the two orthogonal RF fields. For notational 
simplicity, the above cases assume that all pulses are referenced to zero phase. The effect of including non-zero 
phases is simply the addition of the phases of the o> pulse and the ^ pulse to the final phase of the signal. 

Consider serial irradiation, where irradiation at the two frequencies occurs at different times. Using the above 
operators and starting from thermal equilibrium, it is found that a pulse of length r/ at <» followed by a pulse of 
length t p " at co 0 with a delay of rbetween the pulses results in an expectation value oscillating at co t given by 

<U0> = (iV„ 0 -iVi)sin(^ (3) 



where t is the time after the end of the second pulse. N° and N° are the thermal populations of the eigenstates 
i-) and |0> , so that the amplitude of the signal depends on the initial difference in the populations connected by 
the first transition excited. Note that for a single crystal the maximum signal occurs when the first pulse induces a 
nutation angle of Q, ,/ = n/2 and the second pulse has Q 10 = or for a crystal oriented for the most efficient 
excitation y N B u = n/2 and y N B l0 t p * = n. For a powder, the observed signal is proportional to the average over 
all possible crystal directions of the nuclear spin angular momentum projected along the axis of the detection coil. 
Usually no signal at <y + is observed using receiver coils oriented in the x>- y> plane. With a detection coil along the 
z'- axis the signal is proportional to 

mrtAy'l (4) 



2tt n 2.-r 

S t . oc jdajdpsinpjdy 



Performing the powder average numerically, it is found that the maximum attainable signal occurs when y N 5, r/ = 
2.13 rad and y„ B w r> 4.26 rad rather than n/2 and rc ( the aligned single crystal results) because randomly 
oriented crystals will experience RF pulses reduced by the angular factors in Q,.and Q l0 (see Eq. 2). This 
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behavior mimics that seen in single frequency NQR where the maximum signal for a powder occurs at a nutation 
angle of 2.08 rad, a third longer than the nutation angle needed (ir/2) for a properly oriented single crystal. For a 

powder, the three-frequency maximum signal size is ffi'^ -fj^ of the maximum signal of a single- 
frequency NQR experiment at <y + . 

Similarly, simultaneous irradiation of the sample at to. and at co, results in an expectation value oscillating at a>. 
such that 

(/,(/))= sin 2£(l-cos#)sin(aM) 

x[(tf 0 °-N+)(cos 2 S + sin 2 $cos*)-(Ni-Nl)(l + cost?)} 

From examination of the geometrical terms, a maximum signal for a single crystal occurs when £=*/8 and & =k, or 
for a crystal oriented for the most efficient excitation B u IB W = tan(*/8) and Yl fl w 1, = 2tt cos(ic/8). Using (7,(0) of 
Eqn. 5, numerically integrate to find the powder-averaged signal in the detection coil oriented along z'. Again no 
signal at a>. is observed using receiver coils oriented in the / plane. Although the signal is dependent in a 
complicated manner on the thermal populations, and therefore tj, the parameters which give a maximum signal 
depend only slightly on (5,. /B l0 =tan(*/8) and y N B l0 t~l& rad, approximately a third longer than 2k cos(*/8)). 
The maximum signal ranges from 67% to 53% (0<77<1) of the corresponding maximum signal for a single- 
frequency experiment at co + . 

In practice, there is a small distribution of quadrupole interactions within the sample so that H^Hq 0 + A# c . A 
second simultaneous pulse applied at time rafter an initial simultaneous pulse can undo the dephasing caused by 
this distribution so that a spin-echo is formed at a time t = r after the end of the second pulse (see FIG. 5). 
Assuming AH Q is small enough that it has a negligible effect on the time evolution during applied pulses (i.e., AH a 
« H u , H l0 ), it is found | TO > = e ^^e^'' n e'' |TO)> , where K*=(cosf 

7,+sin^ Iy) r corresponds to the first simultaneous pulse of length t p " and ^=(cos^*/ I +sin^/ > ) tf b to the second 
pulse of length t„ b . The refocused signal at <y + is then 

(/, (0) = - sin 2^ (1 - cos & ) sin(fi> + (< + 1\ + r+ 1 ; ) + Aco + (t -t) + 2A<p) (g) 
x[(W 0 °-iV:)(cos 2 ^+sin 2 | fl cos^)-(^!-iV:)(l + cost?' , )]xsin 2 ^cos 2 <r(l-cosi? t ) r 

where Aa>. describes the distribution in co + and A<* « A£+A&. (A£ is the phase difference between the first and 
second pulse of <o. and A& the phase difference for the a* pulses.) For a single crystal, the distribution due to AH Q 
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is completely refocused and the largest echo occurs when <r=ir/8, 6=it/4, and & a =& b =n (or for a crystal oriented 
for the most efficient excitation Y^ Xf f p °=2rc cos(tt/8), ^ 10 f p *=2n cos(rc/4)). Through numerical integration, it is 
found that for a powder the echo is a maximum when B n l JB a l0 *t<in(n/$), jB*,/B* l0 etan(ic/4) f fiv^io t£*lA rad, and 
r*Pvtf*M rad - The maximum signal ranges from 48% to 41% (0<^<1) of the corresponding maximum signal 
for a single-frequency resonant experiment at o) +9 or approximately 75% of the signal following the first pulse is 
refocused for a powder sample. When deriving Eqn. (6) there are several terms not included that depend on the 
precise nature of the quadrupole field distribution and may give rise to echoes at times other than at t = r, as 
discussed by Grechishkin. 



Disclosure of Invention 



The three eigenstates for spin-1 nuclei under a quadrupolar Hamiltonian and the corresponding three possible 
transition frequencies are shown in FIG. 1. In conventional NQR, only one of those transitions is irradiated. The 
resulting nuclear magnetization oscillates at the frequency of the irradiated transition and is observed. Three- 
frequency NQR involves excitation of at least two transitions that causes an observed signal at the third transition 
frequency. Two transitions are irradiated to create an oscillating signal from the third. Similar techniques have 
been used in NMR to create multiple-quantum coherence, but there the excited transition is forbidden and not 
directly observable at its own frequency. 

The ability to detect a signal at a frequency different from the irradiation frequency gives three-frequency NQR 
some important potential advantages over single-frequency NQR. Interfering signals from resonant acoustic 
ringing (magnetostriction) of certain metals at the irradiation frequency could be eliminated. In addition the 
receiver dead-time, the time after sample irradiation in which the receiver's response to the transmitted pulses 
obscures the initial NQR signal, should be reduced since no RF energy is applied at the frequency of the receiver. 
This reduction of dead-time would be especially important for those materials in which the free induction time (7 2 *) 
is short. 



It is an object of the present invention to provide a multiple-frequency technique of nuclear quadrupole resonance 

A further object of the present invention is to provide a system for detecting explosives and narcotics by nuclear 
quadrupole resonance. 

Additional objects and advantages of the invention will be set forth in the description which follows, and, in part, 
will be obvious from the description, or may be learned by practice of the invention. 

Objects of the present invention are achieved by using two frequencies to excite the third transition, and then to 
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detect the signal from this third transition. This avoids any interference from acoustic ringing since the RF 
magnetic field is not applied at the frequency that is detected. The method disclosed here involves the direct 
observation of a NQR transition near a frequency that has not been used for irradiation. 

Brief Description of Drawings 

The present invention will become more fully understood from the detailed description given herein and the 
accompanying drawings which are given by way of illustration only, and thus are not limitative of the present 
invention and wherein: 

FIG. 1 illustrates energy levels determined by a quadrupolar Hamiltonian of nuclear spin 1=1. 
FIG. 2 illustrates a block diagram of the NQR system for the present invention. 
FIG. 3 illustrates a schematic of the three-frequency NQR system. 

FIG. 4 illustrates the dependence of the three-frequency NQR signal of a powder of sodium nitrite at co+ on the 
nutation angle, y^ B 10 r p , for the serial pulse sequence shown in the inset (r = 0). 

FIG. 5 illustrates the dependence of the three-frequency NQR signal at eo + on the nutation angle y N B l0 1 ? for the 
simultaneous pulse sequence shown in the inset 

FIG. 6 illustrates the observed three-frequency spin-echo signal (symbols) as a function of the nutation angle y N 
&*\o t p * for the parameters given. 

FIGS. 7 A & 7B illustrate the problem of acoustic ringing for single-frequency irradiation. 
FIGS. 8A & 8B illustrate the problem of acoustic ringing for simultaneous irradiation. 
FIGS 9A & 9B illustrate the elimination of acoustic ringing for serial irradiation. 

FIGS. 10A & 10B illustrate steady-state pulse sequences for simultaneous (FIG. 10A) and sequential (FIG. 10B) 
pulse irradiation. The sequences are repeated for a total of n signal acquisitions. 

FIGS. 1 1A & 11B illustrate relative SNR as functions of excitation pulse widths for simultaneous (FIG. 1 1A) 
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sequential (FIG. 11B) pulse irradiation where are the pulse widths at 5.19 MHz and are the pulse widths at 
1.78 MHz. The sample is RDX. 



Modes for Carrying Out the Invention 

Reference will now be made in detail to the present preferred embodiments of the present invention, examples of 
which are illustrated in the accompanying drawings, wherein like reference numerals refer to like elements throughout. 



ile resonance 



This invention provides a means for exciting and detecting three-frequency pulsed nuclear quadrupole 
(NQR). One of the main applications is in explosives and contraband detection by NQR. For example, the three- 
frequency technique helps to reduce the receiver dead time in an NQR measurement, which can be particularly 
important if a very high-Q receiver coil is used. In addition, the three-frequency technique provides an even more 
material-specific signature than a single-frequency NQR technique. In general, the three-frequency technique, as 
an extension of traditional NQR, is also useful in some circumstances to better study the behavior and properties of 
certain materials. It is also noted that the three-frequency technique with a known and appropriate material present 
can be used as a simple logic gate since the signal at the third frequency is only observed if both of the two 
excitation frequencies are present 

For a spin-1 nucleus such as ,4 N, there are three distinct transition frequencies as shown in FIG. 1. It is important 
to note that in order to use a multiple-frequency technique, the nucleus to be examined should have asymmetry 
parameter n not equal to 0 or 1, for which the three energy levels are distinct. Fortunately, the I4 N found in many 
explosives has a large asymmetry parameter. As an example, for RDX, i\ = 0.62 and the NQR resonant 
frequencies are v t =5.2 MHz, u. = 3.4 MHz, and u 0 = 1.8 MHz. 

A major aspect here is the use of simultaneous excitation at two frequencies to create the desired signal at a third 
frequency, and also the method and techniques for the creation of spin echoes and other signals analogous to those 
obtained with standard multipulse single-frequency NQR methods. Echoes, for example, are highly advantageous for 
increasing the signal-to-noise ratio per unit time and thereby improving the detection. Multiple echos, such as in the 
spin-lock spin echo (SLSE) sequence, are a common way to significantly decrease the acquisition times for weak 
signals. 



Separating the excitation frequencies from the detection frequency allows the use of low Q transmitter coils to prod! 
large bandwidth pulses while using a high Q detector coil. This will be an advantage in super-g detection of broad 



uce 



lines. 
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FIG. 2 shows a block diagram for the NQR detection system for a preferred embodiment of the present invention. The 
RF source (1), the pulse programmer and RF gate (1) and an RF power amplifier (1) are provided to generate a train 
of RF pulses having a frequency corresponding to one of the three NQR transition frequencies of the type of explosive 
(e.g., all RDX-based explosives) or narcotic desired to be detected. The coupling network (1) conveys this train of RF 
pulses to the irradiating and detecting means (typically a coil). Similarly, the RF source (2), the pulse programmer and 
RF gate (2) and RF power amplifier (2) are provided to generate a train of RF pulses having a frequency corresponding 
to a second frequency of the three NQR transition frequencies. The coupling network (2) conveys this train of RF 
pulses to the irradiating and detecting means (typically a second coil). Likewise, the RF source (3), the pulse 
programmer and RF gate (3) and RF power amplifier (3) are provided to generate a train of RF pulses having a 
frequency corresponding to the third frequency of the three NQR transition frequencies. A coupling network (3) 
conveys this train of RF pulses to a third irradiating and detecting means (typically a third coil). The coupling 
networks ((1), (2), and (3)) also conduct the signal to thereceiver/RF detectors ((1), (2), and (3)) from the coils while 
a specimen is irradiated with the train of RF pulses from all three sources. A central processing unit (CPU) controls 
the RF sources and the pulse programmers and RF gates. The CPU also processes the data and compares the NQR 
signal with a predetermined threshold value. When the predetermined threshold value is exceeded, an optional alarm 
is activated in response to the comparison by the CPU. The receiver/RF detectors, the RF power amplifiers, the pulse 
programmers and RF gates, the RF sources, and the CPU and the alarm may be contained in a console with only the 
coil and the coupling networks being outside of the console. The embodiment shown is one of the most general 
possible applications of three-frequency NQR. 

Although in this embodiment the same means is used for both irradiating the sample with the excitation radiation and 
detecting the NQR signal (e.g., one coil is used for both functions), this is not a requirement of the invention. Separate 
irradiation and detection means (e.g., separate irradiation and detection coils) may be employed if desired. For 
simplicity, common transmitter and receiver coils are used for each frequency in the three-frequency NQR, though 
separate sets of coils could be employed. 

The free induction decay curve is a time domain curve. By performing a Fourier transform on this curve, a frequency 
domain NQR spectrum of the target species may be obtained. 

Although a few preferred embodiments of the present invention have been shown and described, it would be 
appreciated by those skilled in the art that changes may be made in these embodiments without departing from the 
principles and spirit of the invention, the scope of which is defined in the claims and their equivalents. 

The technique utilized according to the present invention is pure nuclear quadrupolar resonance as taught in the 
previously mentioned Buess et al. patent Three frequency excitation and detection may be performed by combinations 
of any means known in the art, for example, a surface coil, such as a meanderline coil or a more conventional "volume" 
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coil such asacylindricalor rectangular solenoid, a tomid, or a Heh^oltz coil. Pie NQR is typically perfonned in zero 
magnetic field: no magnet is required. 

Obviously, many modifications and variations of the present invention are possible in light of the above teachings It 
is therefore.* be understood that, within the scope of the claims, the invention may be practiced otherwise than as" 
specifically described. 



Examples 

The experiment was carried out using "N nuclei in a powder sample of sodium nitrite, NaN0 2 , at room temperature 
Grven sodium nitrite's asymmetry parameter of rj = 0.379, the level spacing is unequal and the transition frequencies 
are distinct (co. /2rr = 4.64 MHz, « /2* = 3.60 MHz, and = 1.04 MHz). The 8 gm sample was placed at the 
center of a probe having three mutually orthogonal pairs of coils (see FIG. 3). RF pulses at co 0 were produced by a 
Tecmag spectrometer, suitably amplified, and coupled into a Helmholtz coil oriented along y', RF pulses at co were 
P^^bya^aencygenemtorandcoupledtoaHelmholtzcoUorientedalong^ 

along z' was used to detect the signal at <u + . The strengths of the RF magnetic fields were measured using a small pick- 
up coil; *,s as large as 2 mT were used for this experiment It was shown that indirect mixing of the frequencies 
resulted in only a negligible RF pulse at co+ (J? 1+ < \yT). 

Each set of cods was parallel tuned to the appropriate NQR frequency using variable capacitors. This frequency- 
selective tuning and the coils' mutual orthogonality provide more than 25 dB of isolation between the transmit and 
receiver coils so that RF leakage of *>_ and co, at the a,, receiver probe is significantiy reduced. The RF leakage of co 
and co, reaching the preamplifier is further reduced by analog filters. Therefore the dead-time could be significantly 
reduced. Receiver dead-times as low as 70 were observed, approximately five times smaller than typical dead-times 
seen in single-frequency experiments. 

After filtering, the NQR signal at co + is heterodyned to co 0 by mixing with RF from the co. frequency synthesizer (see 
FIG. 3) so that the random phase introduced by the pulse at co. is removed. The NQR spectrometer then operates 
normally to detect the signal at co 0 , and all modifications are external to the spectrometer. Approximately 120 
measurements were conducted using simple phase cycling to remove undesired transients after the RF pulses. 

HG. 4 shows the peak amplitudes from the Fourier transform of the free induction decay (FID) at co + in response to 
irradiation of the sample first by co., immediately followed by a pulse of the same length at co 0 . The nutation angles 
on the horizontal axis correspond to a constant B l0 and variation in pulse length between 50 and 300 ,is. To first 
approximation, therelaxation occurring during the pulses was taken into account by considering the signal at the same 
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time after the beginning of the second pulse. Also shown in FIG. 4 are the corresponding theoretical predications using 
the independently measured fields B l0 and B u , and the sensitivity of the receiver coil as determined by a single- 
frequency NQR experiment at v + . (There are no adjustable parameters.) Although the data agree fairly well with 
theoretical prediction for smaller nutation angles ( y N B l0 t f < 2k), they deviate for larger nutation angles. 

The amplitudes of the NQR signals at <u + in response to simultaneous irradiation of the sample by RF pulses at and 
co agree well with theoretical predications over the entire range of nutation angles studied. The symbols in FIG. 5 
correspond to signals taken at the same time after the center of the pulse, in order to take into the account relaxation 
occurring during the pulse length (60 to 500 us). By varying the values of the magnetic fields B l0 and B u we fit 
theoretical predictions to the experimental data. The best fits, shown as curves in FIG. 5, correspond to magnetic field 
strengths approximately 6% higher than measured values, well within experimental error bars. 

For a pair of simultaneous pulses separated by time r, we observed a spin-echo around t - rafter the end of the second 
pulse; see the inset of FIG. 6. Phase cycling eliminated the FID from the second pulse. The signal from t = ronwards 
was Fourier transformed to compare it to a single-frequency experiment. The vertical axis in FIG. 6, was normalized 
to the maximum FID signal for a single-frequency experiment at <y v For a given pair of pulses rwas varied (from 5 
ms down to less than 0.5 ms) and the signal amplitude was extrapolated back to r = 0 using a linear fit to the data. A 
linear fit seems to characterize the data well. The extrapolated amplitudes corresponding to r= 0 are plotted as squares 
in FIG. 6 as a function of nutation angle y H B l0 r p " (r„ a ranges from 180 to 420 /u for FIG. 6 data). For this data the 
ratios of S,s for the two sets of pulses were governed by B" ,. / B a 10 = tan(u/8), B" u I B b l0 = tan(ic/4), and 

As can be seen in FIG. 6, the theoretical predictions of the spin-echo amplitude are in reasonable agreement with 
experimental measurements for these optimal parameters. 

The oscillating magnetic field of the RF pulse produces small changes in the physical dimensions of ferromagnetic 
materials (acoustic ringing) and conversely, these physical deformations produce changes of magnetization in the 
material which can then be detected as an artificial signal. The dashed line of FIG. 7 A shows the resulting response 
from a magnetized paper clip to a single-frequency pulse at coJ2tz = 4.64 MHz. For comparison, the response of an 8 g 
sample of sodium nitrite NaN0 2 to the same pulse is shown as a solid line. Following the same peak in frequency space 
(designated by a star FIG. 7B) it was observed that the response of the paper clip to the RF pulse is linear both in the RF field 
strength and in the duration of the pulse, for pulse widths short compared to the time constant of the acoustic ringing, here ca 1 
ms.. 

Forsimultaneousiiradiato 

acoustic ringing signal was observed at a>+ as shown by the dotted line of FIGS. 8A and 8B. However, for serial irradiation 
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of me same paperclip first with u then with no acoustic ringing at <u + was observed as shown in FIGS. 9A and 9B. The 
pulse parameters used in FIGS. 8A, 8B, 9A, 9B are close to the optimum needed for an NQR signal for each sequence (the 
NQR signals are shown as solid lines). The fact that the acoustic ringing at co t appears for simultaneous, not for serial 
irradiation, ™phes mat me deformation ofme paperclip in respoi^e to me os^ 
component The response fiom the paperclip was Imearm me dijration of me simultaneous pulse 

field strengths B t and B, (shown in FIG. 8B). This suggests that the nonlinear component is proportional to the square of the 
total magnetic fields. Therefore the use of serial irradiation of the sample with u and H and detection at a>. looks like a 
promising means to detect a "N NQR signal, without interference from acoustic ringing. Furthermoie.asshowninHG.8A 
and FIG. 9A. the receiver dead time is significantly reduced over that of a single-frequency experiment 

In order to increase SNR, multiple pulse sequences were used where RF pulses are applied in a time scale that is shorter than 
the time required for the spin system to return to equffibrium Apphcation of such steady-state pulse sequence wim me thn^ 
frequencyNQRmethodhasbeendemons^^ 

and sequential irradiation was tested. The pulse sequences are shown in FIGS. 10A & 10B. The RF excitation used 
frequencies cJ2n = 5.19 MHz, (oJ2n = 1.78 and detection was at coJ2n = 3.41 MHz. The excitation and data acquisitions 
were repeated for a total of « times and the resulting signals from the pulses were averaged together. Different combinations on 
non-alternating and alternating RF phases of 0 and n were used for the pulse pair* of fc\ t„\ and with the best on- 

resonance SNR obtained using non-phase alternating pulses. Using excitation RF fields of B u = 0.26 mT and B 10 = 0.28 mT, 
sigr^amptitudeswereob^ 

are shown in FIGS. 1 1A & 11B. Using the optimized pulse widths for the simultaneous and sequential excitation sequences 
statistical measurements "^repetition of 200 expenWus^ 

to have SNR per unit time approximately two times larger than the sequentially excited sequence, and the overall sensitivity is 
comparable to conventional single frequency multiple pulse sequences used for detection of ,4 N NQR signal in RDX. 



12 



WO 03/076952 

NC 83,803 - PCT 
Inventors: Sauer et al. 



PCT/US02/20513 



Claims 

What is claimed is: 

1. An apparatus for examining a specimen by nuclear quadrupole resonance, comprising: 

a device for irradiating a specimen with a first radio frequency pulse along an x-axis and a second radio 
frequency pulse along a distinct y-axis and receiving along a distinct z-axis a signal from said specimen in response 
to said irradiation, wherein said radio frequency pulses and said signal are defined by the nuclear quadrupole 
resonance frequencies of the nucleus being examined. 

2. An apparatus as in claim 1, further comprising: 

(a) a unit for comparing said signal to a predetermined threshold value; and 

(b) an alarm for signaling when said signal exceeds said predetermined value. 

3. An apparatus as in claim 1, wherein said x-axis, y-axis, and z-axis are orthogonal. 

4. An apparatus as in claim 1, wherein said first and second radio frequency pulses are near to a l4 N nuclear 
quadrupole resonance frequency of a predetermined type of explosive or narcotic to be detected. 

5. An apparatus as in claim 4, wherein said predetermined type of explosive or narcotic to be detected is selected 
from the group consisting of RDX-based explosives, HMX, PETN, TNT, ammonium nitrate, potassium nitrate, 
cocaine, and heroin. 

6. An apparatus for examining a specimen by nuclear quadrupole resonance, comprising: 

(a) a first generator for generating a first radio frequency pulse having a first predetermined frequency and a first 
radio frequency magnetic field; 

(b) a first device for irradiating a specimen with said first radio frequency pulse and said first radio frequency 
magnetic field; 

(c) a second generator for generating a second radio frequency pulse having a second predetermined frequency and 
a second radio frequency magnetic field; 

(d) a second device for irradiating said specimen with said second radio frequency pulse and said second radio 
frequency magnetic field; and 
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(e) a detector for detecting a signal from said specimen in response to iiradiating said 



specimen, 



wherein said radio frequency pulses are determined by the nuclear quadrupole resonance frequencies of the nucleus 
be,ng examined and wherein said first device is a first coil about an x-axis for irradiating said specimen about said 
x-ax JS and said second deviceis a second coil about a distinct y-axis for irradiating said specimen about said y-axis 
and sard detector is a third coil about a distinct z-axis for receiving said signal about said z-axis. 

7. An apparatus as in claim 6, wherein said radio frequency magnetic fields are orthogonal. 

8. An apparatus for examining a specimen by nuclear quadrupole resonance, comprising: 
(a) a first coil about a first axis; 



(b) a second coil about a second axis; 

(c) a third coi] about a third axis; 



(d) a first power source connected to said first coil for generating a fust radio frequency pulse having a first 
predetermined frequency and a first radio frequency magnetic field; 

(e) a seeond power source connected to said second coil for generating a second radio frequency pulse having a 
second predetermined frequency and a second radio frequency magnetic field; and 

(f) a first detector connected to said third coil; 

wherein said first and second predetermined frequencies are defined by the nuclear quadrupole resonance 
frequencies of the nucleus to be examined. 

9. An apparatus as in claim 8, wherein said coils are orthogonal. 

10. An apparatus as in claim 8, further comprising: 



(a) a fourth coil about same axis as said first coil; 

(b) a fifth coil about same axis as said second coil; 

(c) a sixth coil about same axis as said third coil; 
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(d) a second detector connected to said fourth coil; and 

(e) a third detector connected to said fifth coil. 

1 1. An apparatus as in claim 10, wherein said axes are orthogonal. 

12. A method for detecting a class of target species containing quadrupolar nuclei in a specimen by nuclear 
quadrupole resonance, comprising: 

(a) generating a first radio frequency pulse having a first predetermined frequency; 

(b) irradiating said specimen with said first radio frequency pulse; 

(c) generating a second radio frequency pulse having a second predetermined frequency; 

(d) irradiating said specimen with said second radio frequency pulse; and 

(e) detecting a nuclear quadrupole resonance signal at a third frequency in response to irradiating said specimen; 

wherein said pulses and nuclear quadrupole resonance signal are defined by the nuclear quadrupole resonance 
frequencies of the nucleus being examined. 

13. A method as in claim 12, wherein said first and second radio frequency pulses are irradiated simultaneously. 

14. A method as in claim 12, wherein said first and second radio frequency pulses are irradiated sequentially. 

15. A method as in claim 12, wherein said first predetermined frequency is set to equal said third frequency and 
a result said third frequency now equals what was the value of said first predetermined frequency. 

16. A method as in claim 12, wherein said first and second radio frequency pulses are near to a 14 N nuclear 
quadrupole resonance frequency of a predetermined type of explosive or narcotic to be detected. 

17. A method as in claim 16, wherein said predetermined type of explosive or narcotic to be detected comprises 
RDX-based explosives. 

18. A method as in claim 15, wherein said first predetermined frequency is irradiated simultaneously as said 
second predetermined frequency. 
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19. A method as in claim 15, wherein said first predetermined frequency is irradiated sequentially with said second 
predetermined frequency. 

20. A method as in claim 12, wherein said first predetermined frequency is irradiated both sequentially and 
simultaneously with said second predetermined frequency. 
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